Thermal drag revisited: Boltzmann versus Kubo 
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The effect of mutual drag between phonons and spin excitations on the thermal conductivity of 
a quantum spin system is discussed. We derive general expression for the drag component of the 
thermal current using Boltzmann equation as well as Kubo linear-response formalism to leading order 
in the spin-phonon coupling. We demonstrate that aside from higher-order corrections which appear 
in the Kubo formalism both approaches yield identical result for the drag thermal conductivity. We 
discuss the range of applicability of our result and provide a generalization of our consideration 
to the cases of fermionic excitations and to anomalous forms of boson-phonon coupling. Several 
asymptotic regimes of our findings relevant to realistic situations are highlighted. 
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I. INTRODUCTION 

Transport phenomena form a prominent group of prob- 
lems in condensed matter physics. They provide a unique 
information on the excitations and their interactions not 
accessible by other methods**** Recently, thermal trans- 
port by spin excitations in low-dimensional quantum 
magnets has received significant attention due to very 
large heat conductivities found in a number of materi- 
als, for reviews see Refs. H, 0- One may speculate that 
thermal conductivity could be used to probe elementary 
excitations in quantum magnets in a fashion analogous 
to the use of electrical conductivity in metals. 

By now it is well established that integrable one- 
dimensional quantum magnets allow for infinite heat 
conductivity^— Experimentally, however, many spin sys- 
tems rather remote from integrability also demonstrate 
large heat conductivities^ Understanding the role of 
coupling of the spin degrees of freedom to an environ- 
ment, such as phonons and impurities, could be essen- 
tial in this context. Phonons are ubiquitous heat carriers 
along with spin excitations in all quantum magnets. Usu- 
ally, interaction between spins and phonons is discussed 
in the context of dissipation of their respective currents. 
Significant progress has been made herej£~— yet, many 
questions remain open. 

In this work we focus on one such question which is 
rarely addressed: the off-diagonal effect of the flow of one 
of the excitations facilitating the flow of the other 
It is referred to as "spin-phonon drag" , in analogy with 
electron-phonon drag discussed in the thermoelectric 
phenomena in metals and semiconductors^—. 

The second question we address in this work is the rela- 
tion between two distinct theoretical approaches to trans- 
port in a generic coupled two-component system, namely 
the quasi-classical Boltzmann transport theory and the 
Kubo linear-response formalism. Such relations, while of 
fundamental importance, remain unclear between many 
techniques^— devised in the past. For selected prob- 



lems and techniques such correspondence has been estab- 
lished rather firmly ; 15 ! 31 but, to the best of our knowl- 
edge, the comparison discussed in this work has not been 
performed previously. 

Historically, the term phonon drag appears in two 
rather distinct contexts. The first connotation is the neg- 
ative effect of phonons on the electrical conductivity by 
slowing down electrons via processes that are different 
from the direct scattering effectsi 17 i 18 The second, also re- 
ferred to as the Gurevich effect,— ^ is responsible for the 
dramatic deviation of the thermopower in many materi- 
als from the predictions of the "electron-only" theory^ 
In its idealized version^^^ the phonon drag results in a 
substantial heat flow due to adjustment of the momen- 
tum distribution of phonons to that of the electrons, the 
latter being displaced by the electric field. 

In this work, the "thermal-only" analog of the Gure- 
vich effect for a generic spin-phonon problem is consid- 
ered. In this case, only the thermal current is of interest. 
In the standard electron-phonon problem, the thermal- 
only drag is discarded traditionally. This is because the 
thermal conductivity by phonons in metals can usually 
be neglected due to strong scattering of phonons and the 
very large ratio of the Fermi to sound velocity^ How- 
ever, this is not the case in several magnetic insulators 
of current interesl£i£i22r— where the magnetic and lat- 
tice heat conductivity can be of the same order. There- 
fore, the drag between spin excitations and phonons can 
be an important phenomenon. We also note in passing, 
that in contrast the previous electron-phonon problems 
the dimensionality of the spin and the phonon system 
in magnetic insulators can be different. Spin excitations 
may be confined to chains, ladders, and planes. Thus, 
the focus of our study is on the general problem of a 
two-component system and the drag effect in the ther- 
mal conductivity. 

While the spin-phonon problem is our main mo- 
tivation, we consider a generic model of bosonic 
quasiparticle-like excitations, e.g. magnons, coupled to 
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phonons. We derive the contribution of the phonon drag 
to the thermal conductivity in the lowest order of the 
boson-phonon coupling using Kubo and Boltzmann for- 
malisms. We demonstrate that both approaches yield 
identical results for the drag component of the thermal 
conductivity, thus establishing a direct correspondence 
between these methods. We note that despite a signifi- 
cant body of work on thermoelectric phenomena, to the 
best of our knowledge, such results have not been dis- 
cussed before within these two approaches. 

While most of the work is devoted to bosonic spin exci- 
tations, we have also generalized our consideration to the 
case of fcrmionic excitations, as well as to the case of par- 
ticle non-conserving boson-phonon interactions. The lat- 
ter are as common as the "normal" , particle-conserving 
ones and occur, e.g., in the phases with broken symme- 
try. We also note that using our expressions for the drag 
thermal conductivity we reproduce some of the results of 
the recent work,— which considers spin-phonon drag in a 
particular class of quantum magnets using the memory- 
function approach. 

The paper is organized as follows. In Sec. |TI]we intro- 
duce the model. Section UTTl outlines the derivation of the 
drag conductivity from the Boltzmann equation. In Sec- 
tion IIVI we detail the Kubo diagrams for drag and elabo- 
rate on one approach to their evaluation. Sec. fVl extends 
the consideration of the drag thermal conductivity onto 
fcrmionic excitations as well as onto anomalous (particle- 
non-conserving) boson-phonon coupling. In Sec. I VII we 
provide a qualitative discussion of thermal drag for sev- 
eral representative cases and asymptotic regimes. We 
conclude with Sec. IVIII Several Appendices are pro- 
vided. In Appendix [X] we detail more technical points of 
the Boltzmann approach. Appendix [Bl is devoted to the 
evaluation of the Kubo diagrams for the drag using an 
alternative approach. Finally, in Appendix [Cl we discuss 
some corrections beyond the Boltzmann results from the 
Kubo diagrams. 



choice of the statistics would be that of fermions. 

The spin-phonon coupling occurs because lattice dis- 
placements cause changes in the spin interactions and 
anisotropics Thus, the simplest, yet very general form 
of the spin-phonon coupling is linear in the lattice dis- 
placement and quadratic in the spin operators. Af- 
ter mapping spins onto bosonic quasiparticlcs, the re- 
sultant lowest-order boson-phonon coupling will gener- 
ally contain the "normal" part, which conserves the bo- 
son number^ and the off-diagonal, "anomalous" bosonic 
terms^ Since the subsequent drag conductivity deriva- 
tion is conceptually identical for the normal and anoma- 
lous forms of interactions, we will postpone consideration 
of the latter until Sec. [V] and will treat it in less detail. 
Here and in the next two Sections, we will concentrate 
on the normal form of the coupling. 

Altogether, the Hamiltonian implied for our subse- 
quent consideration is 

H = -ff b + H p h + -ffb-ph , (1) 

H h = ^k&lX , H ph = ^ w P a P a p . ( 2 ) 
k P 

ffb-ph = £ v^f_ p bt-M a l + a - P ) . ( 3 ) 
k, P 

where and tip ^ are boson and phonon operators, and 

^-p;k- P ,k = (^p b ;k P k- P )* due t0 hermiticity of i/ b -ph, 
and we do not specify interaction terms that result in 
the relaxation rates of bosons and phonons. Note that, 
aside from a more general momentum dependence, the 
boson-phonon interaction in @ is the same as the one in 
the electron-phonon coupling case. 

III. BOLTZMANN APPROACH 
A. Thermal drag conductivity 



II. MODEL 

Spin systems can yield a wide variety of excitations, 
such as spinons, magnons, triplet excitations, etc., de- 
pending on the dimensionality, spin value, and geome- 
try of the structural arrangement. Because the focus 
of this work is on the drag between spin excitations and 
phonons, we assume that both of them are describablc by 
wcll-dchncd quasiparticlcs. For the spin excitations we 
assume energy dispersion and some phenomenological 
intrinsic or extrinsic transport relaxation rate. This rate 
may include scattering caused by impurities, spin-spin 
interaction, and other degrees of freedom. The corre- 
sponding energy for phonons will be denoted by oj p and 
they are also assumed to have a finite relaxation rate. 
In most of the paper, we assume the statistics of the 
spin excitations to be that of bosons. Sec. [V] outlines 
the changes in the drag conductivity which results if the 



Let us denote boson and phonon distribution functions 
as /k and n p , respectively. In Boltzmann's approach the 
total heat current is the sum of the currents from each of 
the particle species: 

jtot = £ v k £k£/k + £ u p o;p(5rip, (4) 

k P 

where Vk = cfek/<9k and u p = dui p /dp are the velocities, 
the chemical potential is set to zero, and <5/k and Sn p are 
the non-equilibrium parts of the distribution functions of 
the bosons and phonons, respectively. 

The distribution functions are determined from the 
Boltzmann equations 

^ = Stl[f,n], ^ = S**[n,/], (5) 

where ^ b ( ph ) are the collision integrals which include all 
possible scatterings for bosons (phonons). These Boltz- 
mann equations are coupled because the boson-phonon 
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interaction in (J3j) yields terms in the collision integrals 
which depend on both / k and n p . 

Assuming the system to be in a steady state under 
a small uniform thermal gradient, we may linearize the 
Boltzmann equations in Sf and Sn to find 



£ k df° _ 6f k ^ 

~T ~de ~~ r£ ^ 



5n r 



ph— >b ' 
P T p,k 



t*, dn° 



b-s-ph ' 
k T k.p 



(6) 



(7) 



where /£ = [e £k / T - l]" 1 and n p = [e"p/ T - l]" 1 are 
the equilibrium distribution functions and the collision 
integrals in the right-hand sides are expanded in Sf and 
Sn and are considered in the relaxation-time approxima- 
tion. The first terms on the right-hand sides of (|6]) and 
([7J arc the usual diffusion terms, with and l/r ph 

being the transport relaxation rates of the bosons and 
phonons due to all possible relaxation mechanisms, as 
discussed in Sec. [TTJ The second term on the right-hand 
side of the Boltzmann equation ([6]) for bosons with the 
momentum k is from the expansion of the collision inte- 
gral St^[f, n] in the phonon non-equilibrium distribution 
Sn p . Because of that it contains an integral over the 
phonon momentum. The same is true for the phonon 
Boltzmann equation ([7]). These latter terms arise solely 
due to the boson-phonon coupling ([3J and are due to the 
non-equilibrium components of the particles of opposite 
species. Therefore, it is natural to identify them with 
the drag from one species of excitations onto the other. 
Below we are going to explicate the relation of 1 /Tp 1 ^ 13 

and 1/Tk^ ph with V^ ) .^^ c _ p through the boson-phonon 
collision integral, but at this stage we simply use them 
as a short-hand notations for the "drag rates" .— 

In general, Eqs. (|6]) and flTJ) reduce to integral equa- 
tions for Sf and Sn. However, we assume that the drag 
terms in (|6]) and ([7]) arc small compared to the diffusion 
contribution, i.e., the drag rates l/r ph ^ b and l/r h ^ ph 
are small compared to the intrinsic boson and phonon 
rates l/r b and l/r ph . This is equivalent to treating 
the boson-phonon coupling V h ~ ph as a perturbation. In 
turn, one can solve and ([7]) iteratively by using the 
diffusion-only components in the integrals containing Sn 
and Sf. This corresponds to neglecting the terms of or- 
der 1 1/ b_ P h | 4 and higher. Within this approximation, the 
total current in ^ is 

jtot = jb +jph + jd,ph->b +jd,b->ph, (8) 

where jb and j p h are the usual, "diagonal" terms, and 
jd,ph->b and jd,b^ph are the currents due to the drag of 
phonons on bosons and vice versa. The total drag current 
can be written as: 

if - i a 4- i a 

-/drag .vd.ph— >h ~ Jd.b— fph 



dn p vi 



k,p 



a P 



duj T ph .^ h de -r-. b ~' ph 



'P,k 



'k, P 



VbT, 



(9) 



where a and /3 are vector components. Choosing the 
temperature gradient in the x-dircction and assuming the 
conductivity tensor to be diagonal, we obtain the drag 
thermal conductivity: 



^drag 



k, P 



dn p 



1 



i 



duj T ph .^ h de r b ^ ph 



'p.k 



'k,p 



(10) 



The above analysis thus far has been independent of the 
microscopic form of the boson-phonon coupling. 



B. Microscopic consideration. 



The "drag rates" l/Tp 1 }^ and 1/T k ^ ph are obtained 



k. P 



by taking variations of / and n in the corresponding func- 
tionals St[f,n]. We now detail the derivation of one of 
them. The collision integral for phonons scattered off 
bosons via Eq. ((3]) contains two terms: the first one in- 
creases the number of phonons with momentum p, the 
second one reduces it. They can be grouped together as: 

Stf[n, /] = 2ttE |Vp b :k p k h - p | 2 • S(s k e k _ p - u p ) 

k 

x [/k(/k- P + rip + 1) - /k- P n p ] (11) 

This is the complete expression of the phonon col- 
lision integral due to phonon-boson interaction ([3]). 
The subsequent linearization of (|lll) uses the condition 
St ph [f°, n°] = 0. Writing f = f° + Sf and n = n° + Sn 
and neglecting terms of order SfSn and (Sf) 2 yields the 
first and second terms in the right-hand side of Eq. (J7J . 
The first one is the diffusion term while the second one 
is the drag term. In the latter, for the terms contain- 
ing (5/k-p, we shift summation over k — > k + p so that 
Sfk-p —> Sf k . After these manipulations, the drag rate 
of bosons on phonons is given by 



St ph [n°J° + Sf}^-J2 



Sfu 



b-ph ' 

k 'k,p 



1 



b-s-ph 



= -2tt V 



'k, P 



rb-ph 



p;k,k— p I 



\fL P + n° p + l] (12) 



x <5(e k - fk-p 



K 



b-ph 
p;k+p,k I 



k+p 



^kH 



- ! P ) 



P 

= Ik 



the second 

b-ph 1 2 



Note that under p - 

changes to \V*~^_ p ^ 
velocity changes its sign, u p —r u_ p , 
tcr ([3]). Using this symmetry we obtain a compact form 
for the component of the thermal conductivity due to the 



p;k,k— p I 



term in ifTZj) 

and the phonon 
sec discussion af- 
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drag of bosons on phonons: 



IV. KUBO APPROACH 



K 



b->ph _ 
drag 



2 7; 



b-ph 
p;k,k— p | 



k. P 



-£r [/k- P + n l + !] %k - £k- P - Wp) (13) 
^[/k°- p -n p ]%k- £ k- P +- F 



The derivation of the drag rate 1 / Tp and of the drag 
conductivity of phonon on bosons follows similar reason- 
ing and is presented in Appendix[5]l. After some algebra 
one arrives at the following statement: 



b— >ph ph— s-b 

K drag = K drag ■ 



(14) 



This relation bears a simple and general physical mean- 
ing: within linear response, the non-equilibrium compo- 
nent of one species of particles causes the same drag on 
the other species as the non-equilibrium component of 
the other causes on the first one. That is, phonons drag 
bosons the same as bosons drag phonons. Thus, the total 
contribution to the conductivity is simply twice the con- 
tribution in Eq. (JT3J. In addition to the algebra above, to 
obtain (|14p we have used the following identities between 
the combinations of the bosonic distribution functions 
and their derivatives, 



dn°„ 



f oi — d/k r to 

'kj - ~0^~IVk-p 



1] 



[fl 



f° 

duo LJK - p de 



k [/.° 



k P 



(15) 

P 

(16) 



which can be obtained with the help of Stp^[f°, n°] = 0. 

Thus, within the Boltzmann formalism, the total drag 
thermal conductivity, to leading order in the boson- 
phonon coupling, is given by 



4 "£ 

k, P 



8fk [ rO 

Uk-p 



de 



p;k,k— p | 

f ri° + 1] 6 (e k - e k _ p - Wp) (17) 
n p ](5(e k - e k _ p + uj f 



This expression is the main result of this work. We 
would like to note, that within the approximations dis- 
cussed above this result is valid for any type of scattering, 
impurity, boundary or Umklapp, all being implicitly in- 
corporated in the transport relaxation times of phonons 
and bosons. This expression also contains normal as well 
as the Umklapp boson-phonon scattering. That is, the 
quasimomcnta in (|17[) are defined up to the reciprocal 
lattice vectors and the summation over the latter is as- 
sumed as usuali^ 



A different theoretical approach to transport, alter- 
native to the Boltzmann equation, is the Kubo linear- 
response formalism. The great advantage of this ap- 
proach is its conceptual clarity with regard to the def- 
inition of the drag thermal conductivity. It is also very 
effective in classifying terms by their respective order in 
the coupling constant as it contains them explicitly. 

In Kubo's approach the uniform part of the thermal 
conductivity is obtained by taking the DC-limit of the 
imaginary part of the dynamical heat-current suscepti- 
bility X ^^jl5 



lim 

ix>->0 oj + iO 



Im [x^ (0,oj + K 



(18) 



where fi and v are the spatial directions and the suscep- 
tibility is a sum of diagonal and off-diagonal terms, 



= E 



(lis 



(19) 



with the components 

X&(q,w + i0)=i/ (L7mi- q7 -]>e i( " +i0)t *, (20) 



which contain the heat-currents j^. In this study, the 
long wavelength limit of the thermal current of bosons is 
given by 



Jql 



k 



e k w k &. 



and the phonon one by 



Jq2 



(21) 



(22) 



where the energies and velocities were defined previously 
in @ and |QJ. For the remainder of the paper the usual 
limit of q = is implied for the currents, however q is 
kept visible for clarity. Note that Eq. (fTS)) is derived from 
the linear response to the temperature gradient, which 
couples to the total energy density. Therefore, apart from 
the bare heat currents of (|2"Tj) and (|22|) . the interaction 
term in the Hamiltonian ([3]) will also give rise to a con- 
tribution to the thermal current. This current, labeled 
by j3. q , follows from the continuity equation, 



q-j3, q = [H,Hq] -q - (jl, q +j 2 ,q), 



(23) 



where = ^ r e _lq ' r i? r is the Fourier transform of 
a position-dependent Hamiltonian energy density, H = 
H r . However, j3 iq does not constitute a contribution 
to thermal drag and we will not consider the correspond- 
ing terms in this work. 

We would like to emphasize that (|18|) refers only to the 
DC-limit and does not incorporate the Drude weight. 4 
The latter is assumed to be zero henceforth. 
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(A) b-ph (B) b ph 

v p;k,k-p v p;k.k-p 




,b-ph * i #b-ph * 

p;k,k-p » p ;k ,k-p 



To calculate the thermal conductivity one needs to per- 
form the frequency summations in (f2"5j) and ([25)1 . We 
utilize two technical approaches for that: the first uses 
the spectral representation for the Matsubara Green's 
functions and the second one uses integration along the 
branch cuts of the Green's functions.— Below we elab- 
orate on the use of the first one while the branch cut 
integration approach is discussed in Appendix [B] 



FIG. 1: Graphical representation of the lowest order, off- 
diagonal current-current correlations contributing to boson- 
phonon drag conductivity. Solid lines are bosons, wavy lines 
are phonons. e's are the auxiliary frequencies used in the 
spectral representation approach. 



The diagonal (i = j) components in (j2"0)l are the "usual" 
diffusion terms and they do not contribute to the drag. 
Naturally, the drag is given by the off-diagonal current- 
current correlation functions, xi,2 and X2,i- Considering 
the boson-phonon coupling 7 p .^_ p in §3§ as a perturba- 
tion, the lowest-order diagrams contributing to the drag 
are shown in Fig. [TJ These two diagrams are the only 
"drag" diagrams of the order |l/ b ~ ph | 2 that contribute 
to Xi,2- The mirror-reflection of these diagrams with re- 
spect to a vertical line yields equivalent contributions to 
X2.i- This is, again, a graphical way of stating that in 
the linear response the drag of phonons on bosons and 
the one from bosons on phonons are identical. Thus, 
the total drag conductivity is given by twice the value 
of the diagrams in Fig. [T] The choice of the momenta 
in Fig. [T] is made to keep all the momentum-dependent 
functions between the two diagrams, such as energies and 
vertices, the same. Since in the diagram A the phonon 
momentum is — p, its sign is opposite due to the velocity 
in the current vertex. The analytical expression for the 
total drag thermal conductivity given by Fig. Q] and its 
mirror-reflection is: 



Kdr 



lim 



Toj ^ 
k. P 



■£k) 



V 



b-ph 



p;k,k— p I 



Im 11a (k, P, w + iO) - n B (k,p,w + iO) 



(24) 



with 



A. Spectral Representation Approach 

The spectral representation for the Matsubara Green's 
function is: 



A k (e) 



IU}„ 



(27) 



where the shorthand notation J £ = J^° oo de/2ir and the 
following relation of the spectral function to the retarded 
Green's function, Ak(e) = —21mG^(e), are used. 

Below we calculate the contributions to the thermal 
conductivity from the diagram A in Fig. [T] those from 
the diagram B are discussed near the end of the section. 
Using the spectral function representation (|27[) . and as- 
signing auxiliary frequencies according to the diagrams 
in Fig. [T] one can rewrite as: 

iTi(k,p,iu;) = f ^k(ei)^k(£ 2 )A k „p(£ 3 ) (28) 

Jei,...es 

xA p (£ 4 )A p (e 5 ) ■ n j i(iw,ei...s), 

where £i...5 stands for the five frequencies associated with 
each individual line in Fig. [TJA). The frequency summa- 
tion over wi and 0J2 is now accumulated in n^, which is 
given by: 



-1...5 



1 



lLO\ — £1 IUJ\ + VjJ — £2 



iuj\ + iu>2 — £3 — £4 — iio ~ £5 



(29) 



Performing Matsubara frequency summations in (|29[) we 
obtain: 



n A (k, p, iw) = T 2 ^ G^(iLu + iwi)G k (iwi) (25) UA ( iuj > £l -s) 



1 



1 



x G k _ p (MJi + iu} 2 )G p (iuj2 - iw)G p (zw 2 ), 



and 



iuj + £1 — £2 lOJ + £5 — £4 
(n 3 - ni)(n 3 _i - n 5 ) _ (n 3 - ni)(n 3 _i - n 4 ) 
lOJ + £5 + £1 - £3 £4 + £1 - £ 3 

(n 3 - n 2 )(n 3 _ 2 - ra 4 ) (n 3 - n 2 )(?i 3 _ 2 - n 5 ) 



(30) 



n B (k, P , iu) = T 2 G *( IUJ + *wi)Gk(*wi) (26) 

x Gk- P («wi - iu} 2 )G p (iuj2 + iuj)G p (iuj2), 

where G^iu) and G p (iuj) are the boson and phonon 
Green's functions, respectively. 



-iuj + £4 



£3 



£5 + £2 - £3 



where 



n°(£j) are the Bose distribution functions 



with the corresponding energies and nj_j = n . (gj — £j). 

For the uniform, DC thermal conductivity (|2"4"|) we need 
to take the imaginary part and the to — > limit of (|30l) : 
ImHA(it-o n — > w + iO) at u — > 0, which splits naturally 
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into four terms: 



Imn A (w) = ^ImlT^V) 



(31) 



= Im (I) Re (II) Re ([III]) + Re (I) Im (II) Re ([III]) 
+Rc (I) Re (II) Im ([III]) - Im (I) Im (II) Im ([III]) , 

where I, II, and [III] are the first, second, and third fac- 
tors of the product in ([3T)| . respectively, and [III] includes 
all the terms inside the square bracket. In what fol- 
lows, we refer to the four contributions to the conduc- 
tivity coming from the four terms in pip as to Imll^™ , 
m = 1 . . . 4. 



1. Boltzmann terms 

Here we explicitly evaluate the leading-order contri- 
butions to the thermal drag conductivity and show that 
the Kubo approach yields the same answer as the one 
obtained using Boltzmann equation. The discussion of 
the other, subleading non-Boltzmann contributions is de- 
ferred to Appendix [Cl 

We would like to assert that within the spectral repre- 
sentation calculation, the leading contributions are given 

only by Imll^ in ([3T|. The rest of the terms, Imn^, 

Imll^ , and ImII\ , yield results that are subleading 
in the sense of containing higher power of r's, which is 
equivalent to having higher-order terms in V h ~ ph and 
other couplings. 

Consider Imll^ term where all three factors in ([30)) 
contribute their imaginary parts, 

IraH^ (u) = it 3 S(uj +£i — s 2 )5(u + e 5 - e 4 ) 

[(n 3 - m)(n 3 -i - n 5 )5(uj + e 5 + £i - e 3 ) (32) 
- (n 3 - n 2 )(n 3 _ 2 - n 4 ) S(e 4 + e 2 - e 3 - u)] . 



parts of their self-energies, r k and s p , such that ru{s p ) -C 
£ k(<j-'p). They are also related to the relaxation times 
used in Scc. lIIII as ty 



2r£ and s" 1 



2T-P 11 , seeRef.[Tf 



Thus, the spectral functions of bosons and phonons can 
be approximated as Lorcntzians: 



A k (ei) = 



2r k 



(e s - cj p ) 2 + s p 



(35) 
(36) 



Since the spectral functions ([331) an d (|3"6"[) arc strongly 
peaked at £k and w p , the main contributions in the in- 
tegrals in (|3"4"[) are obtained at E\ w £k and £5 « oj p . 
Identifying distribution functions with that of phonons 
and bosons via n(u p ) = n p , n(e-k) = /£, an d "(e'k-p) = 
/k-p> Anally yields, 



Im IlA(k,p,w) 



ui 



;=0 



k p "&7 L/k- P - /kj 



x A k - P (ek + w P ) 



(37) 



- £k-p), 



where in the last line we have approximated the Loren- 
zian with the delta-function and have neglected contribu- 
tions from 11^4 3 ^ terms. As we discuss in Appendix [Cj 
both approximation are of the same order and correspond 
to neglecting terms that are subleading to ([37)) . 

Repeating the same consideration for the diagram B 
in Fig. [JJ gives: 



Im II s (k,p,w) 



u=0 



9 b P h r ,0 _ f oi 
/7rr k T p duJ L/k-p /kj 

x 5{e k - uj p - £k- P ) • (38) 



In the low frequency limit Eq. (|32J) reduces to 



imir^V) 



u=0 



2tt 3 5(s 1 - s 2 )S(e 5 - e 4 ) 



x (n 3 -ni)— 5(e 5 +£i-£ 3 ). (33) 

Substituting this into (|28[) and performing integrations 
with delta-functions, we obtain the leading contribution 
to ImhL^k, p,w), 



Im4 4) (k,p, W ) 



j=0 



A k (£i) 



^p(£5) 



9n B 



-4k- P (ei + £5) ["1+5 - "1] -gj- ■ (34) 

We assume that bosons and phonons are well-defined 
quasiparticlcs with frequency-independent imaginary 



Using the identities for the distribution functions in (|15l) , 
([TB)l and substituting j37| and ((38j) into (JUJ gives the 
Kubo answer for the drag component of the thermal con- 
ductivity 

Kdrag = -f £« £k7 k b ) (^ P fl |v p b ;k tp| 2 



k. P 



[/k°- 



p + n p 



p + l]5(£k-£k-p-W p ) (39) 



[/k-p-np]*(ek-ek-p + w P ) 



One can see that this is identical to the Boltzmann an- 
swer in ([TT| . We discuss in Appendix [Cl that the con- 
tributions from the remaining terms Imn 



= 1...3 

and corrections to ([39| due to the broadening in the spec- 
tral functions are of the order 0(t 3 ) and, therefore, can 
be neglected. 
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One can check the consistency of the drag conductiv- 
ity expression in (fT?]) and with the diagonal terms 
in conductivity by assuming that the leading source 
of the relaxations defining both the spin and phoiion 
transport relaxation times is the spin-phonon coupling 
in <j3j> - Then the diagonal and the drag conductivities 
are all of the same order in the spin-phonon coupling: 
K ph _ K b _ Kdrag K i/|v"b- P h|2_ Thig alsQ demonstrates 

that in the idealized case of free spin excitations coupled 
to dissipationlcss phonons via a weak coupling all con- 
ductivities should be of the same order in that coupling. 

We would like to emphasize again, that the identity 
between Eq. ([Ml) and Eq. ([T7]) is rather remarkable as 
they are derived starting from completely different phys- 
ical formulations. 



V. FERMIONS AND OTHER 

Here we generalize the analysis of this work onto two 
additional cases. First, we assume the same form of the 
coupling to phonons ([3]), but consider fermions instead 
of bosons. This scenario is not only applicable to cases 
where the spin algebra has been mapped onto fermions, 
but is also relevant to the thermal conductivity in met- 
als and semiconductors. The second generalization ex- 
tends our drag consideration on the case of anomalous 
bosonic terms in the boson-phonon interaction. Such 
terms readily exist in the interaction of phonons with 
magnons in the ordered antifcrromagnets, as was dis- 
cussed previously^ They also exist for triplet excitations 
in gapped, dimerized, and other phases. The following 
derivations are based on the Boltzmann formalism only. 



A. fermions 

Coupling of fermionic excitations with phonons is, gen- 
erally, of the same form as given in Eq. §3§ . This is obvi- 
ously the case for the coupling of electrons with phonons, 
and is also true for the XX Z spin chains when spins are 
represented by the Jordan- Wigner fermions. The gen- 
eral expression for the drag current will be still given by 
Eq. @, where £k is replaced by £k = £k — the fermion 
energy relative to the chemical potential, and the "drag 
rates" 1 /r plwf and 1 / r f ^ ph are determined by the corre- 
sponding collision integral involving fermions and bosons, 
with / now representing the fermion occupation number 
and = (cxp[(£k — fi)/T] + being the equilibrium 
Fermi-distribution function. One obvious difference for 
the probabilities is that a fermion with the momentum 
k is created with the probability given by (1 — /k)- The 
derivation for the drag conductivity follows exactly the 
same steps as those for the boson-phonon case considered 
in Section [IIII Useful identities for certain combinations 
of f° and n°, analogous to the ones in (jT5j) and (TTS)) . are 
listed in Appendix \K\2. Taking into consideration the 



above differences we obtain the total drag conductivity 





47T 


^drag — 


~T 


X 






de 



r ph\ lyf-ph I 
p / r p:k,k-p I 



k,p 



-B [1 - /k°- P + S(s k - £k - p - up) (40) 



df° 



To summarize, the phonon drag conductivity for the 
fermionic case (|4"0")) takes the same form as for the bosonic 
case (fT7"|) with two modifications: (i) / k _ p — > ^ik-pi ( u ) 
£k — > Ek = £k — Note that the second change should 
also be made in the case of bosons if the chemical poten- 
tial for them is not zero. 

We note, that the fermionic case of the drag discussed 
here is different from the one traditionally considered 
in the thermoelectric phenomena. As is mentioned in 
Sec. U for the electron-phonon problem in metals, the 
thermal-only drag effect is usually neglected because of 
the dominance of the electronic thermal conductivity over 
the phonon onci This is not the case in many low- 
dimensional quantum magnets &2£2r£5. 

Regarding potentially different outcomes of the drag 
effect for the fermionic systems (|4T)]) compared to the 
bosonic ones (fT?]) [(|39p], we remark that the major differ- 
ence may arise due to the presence of the Fermi surface 
in the former cases. It is known, that the "normal" and 
the Umklapp scatterings contribute with opposite sign 
to the drag conductivity, which is discussed as one of 
the reasons for the suppression of the Gurevich effect in 
metals^ ' 16 i 20 Such an effect of Umklapp can be expected 
to be small at low temperature for the bosonic case be- 
cause all of the heat carriers are at small momenta. For 
the fermionic case, on the other hand, the effect of the 
Umklapp should be present, similarly to the electron- 
phonon case. However, since the fermionic representa- 
tion of spins is restricted to ID, significant differences 
from the traditional 3D electron-phonon consideration 
may also occur. Any quantitative statement on whether 
the drag will be more substantial for magnetic excita- 
tions obeying bosonic or fermionic statistics will depend 
on specific model calculations, which are not the focus of 
this work. 



B. anomalous bosonic terms 

Next we consider drag contributions in the phonon- 
boson case due to anomalous terms of the kind 



^ = E^ b k p k-p ()Wk«P+H.c 
k, P 



(41) 



These describe processes involving creation of two bosons 
from a phonon and generation of a phonon due to the 
annihilation of two bosons. 
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With the details of the algebra provided in Ap- 
pendixes, here we simply state that the approach de- 
scribed in Sec. Hill yields the following result 

k,p 

f)n° 

x a? [l + /k ° + a(£k + £k - p ~ Wp) ■ (42) 

In the case when both the "normal" and "anoma- 
lous" (|1T|) boson-phonon couplings are present, the lead- 
ing contribution to the drag thermal conductivity is the 
sum of the results in (jTT)) and (O . 

VI. QUALITATIVE ESTIMATES 

In this section we provide a qualitative discussion of 
various asymptotic results that can be readily inferred 
from Eq. (|17[) for several representative spin-phonon sys- 
tems with the goal of estimating when drag effects can be 
significant and when they are not. The temperature de- 
pendence of the drag thermal conductivity is determined 
by two factors: scattering lifetimes and the occupation 
numbers of the excitations. 



A. Boundary-limited regime 

First, we would like to consider gaplcss spin excita- 
tions with linear dispersion £k ~ v\k.\, coupled to acous- 
tic 3D phonons, the situation relevant to a wide vari- 
ety of antiferromagnets. For the low impurity concen- 
tration and at low temperatures both phonon and bo- 
son mean-free paths can be expected to be boundary- 
limited, the case well documented for Nd2Cu04^ How- 
ever, the heat carrying excitations will be few in num- 
ber and the drag conductivity has to go to zero at low 
temperatures. A straightforward algebra in (fTT)) yields a 
power-law: K<jrag oc T 7 , with 7 = 2 + D s + m, where D s 
is the dimensionality of the spin system and m depends 
on the long-wavelength k- and p-dependence of the spin- 
phonon coupling V^,.^_ . In the case of D s = 3 (e.g., 
3D magnons) and assuming that the coupling follows the 
standard form V^. k p ^i oc y/pkk', which corresponds to 

m = 3, altogether gives K|jrag K ^ '• This should be com- 
pared with the diagonal thermal conductivities in this 
regime « ph ~ n h oc T 3 . Thus, the drag effect is, gener- 
ally, subleading in the considered regime. 

B. Gapped spin system 

In another specific example let us consider a gapped 
spin system at low enough temperatures so that the occu- 
pation number of spin excitations is exponentially small: 
/k oc e~ A / T , where A is the gap in the spectrum. In 



the case when the relaxation within the spin system is 
only due to a weak coupling to phonons whose relaxation 
rate is dominated by the Umklapp processes, i.e. T ph oc 

e e °/ T , where 0£> is a fraction of the Debye energy^ 
our Eq. (jTTJ) naturally leads to K dr a g oc e (§D_A)/T . This 
result was obtained in Ref. [IH using the memory-matrix 
approach. 

C. High-temperature regime and disorder effects 

Third, we consider the high-temperature limit, for ei- 
ther gapped or gapless spin system, when temperature 
is higher than both the Debye energy and the spin- 
excitation energy scale, T ^> @d,J- Formally this case 
may raise questions regarding the transition into the dis- 
ordered state, however is fully analogous to the textbook 
consideration of the lattice thermal conductivity at T > 
0£>ii In fact, in this region quasiparticles can be con- 
sidered as strongly damped. The rates of the Umklapp 
scattering for spin excitations and phonons are high and 
are proportional to the occupation numbers of a "typi- 
cal" boson or phonon, thus leading to r ph ~ r b oc 1/T. 
The rest of the estimate in Eq. (fP7|) is again straightfor- 
ward, giving Kdrag oc 1/T. This should be compared to 
the diagonal conductivities in this regime, which show 
the same asymptotic behavior^ « ph ~ n h oc 1/T. This 
consideration, combined with the low-temperature one, 
implies that the drag conductivity should go through a 
maximum at intermediate temperatures, similar to the 
diagonal conductivities. 

When the energy scales of the phonon and spin sys- 
tem are well separated, as in the cuprate-based materi- 
als where J ^> Od, another asymptotic regime is pos- 
sible, 0£> <C T <C J. Intuitively, the drag can be ex- 
pected to diminish together with the phonon conductiv- 
ity (« ph oc 1/T) because phonons are sufficiently equili- 
brated by the phonon-phonon scattering. However, the 
T-dependence of the drag also depends on the specifics of 
the relaxation within the spin system. Thus, no definite 
conclusion on the prevalent behavior of the drag conduc- 
tivity in this regime can be drawn without identifying 
such a relaxation. 

The disorder dependence of the drag can also be con- 
sidered using similar qualitative reasoning. If the dis- 
order affects both types of excitations on equal footing, 
so that K ph ~ K h oc l/ni m p, where ni mp is the impurity 
concentration, then the drag conductivity diminishes as 
Kdrag oc l/(ni mp ) 2 . If the disorder can be introduced se- 
lectively in one of the sub-systems without significantly 
affecting the other, as in the case of lattice disorder in the 
ladder cuprate system CagLa5Cu2404i j£ the drag con- 
ductivity will be reduced together with the diagonal con- 
ductivity of the most affected species of excitations. 

Thus, intuitive conditions for maximizing the effect of 
drag arc the simultaneous presence of significant popu- 
lation of spin excitations and phonons with long scatter- 
ing times. Since such conditions also imply large diago- 
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nal contributions of spins and phonons to the heat cur- 
rent, they are typically satisfied for temperatures that 
are low enough in comparison with cither J or 6d 
but are above the boundary-limited regime. Note that 
the optimal regime for the phonon drag in thermoelec- 
tric phenomenon is often quoted as T ~ Op/ 5.— Such 
a regime can be of relevance to the recently reported 
record-breaking thermal conductivity by spin excitation 
in a high-purity ID spin-chain material SrCuC>2, Ref . l35l. 
where a nearly ballistic propagation of spin excitations 
was reported. 

The issue of the separation of the drag component of 
the thermal conductivity from the "diagonal" one may 
require a series of doping experiments in which disorder 
is introduced deliberately to suppress the conductivity 
of one of the species and thus diminishing the drag as 
welliM 



D. qualitative estimate of the drag 

Lastly, we would like to come back to the problem of 
the gapless spin excitations with linear dispersion cou- 
pled to phonons, the problem motivated by the ID spin- 
chain and 2D layered cuprates where the spin excitations 
are fast and the phonons are slow, J ^ <3d ■ Analo- 
gous to similar estimates of the thermoelectric power^ 
and without reference to a specific model, the following 
consideration is not intended to be entirely rigorous, but 
rather is aimed at deriving an upper-limit estimate of the 
thermal Gurevich effect. 

Let us assume that the boson relaxation is due to 
impurities or some other extrinsic or intrinsic mech- 
anism while phonons are dissipationlcss, a considera- 
tion similar to the electron-phonon drag problem^ Such 
a scenario is also potentially relevant to the ID spin- 
chain materials in low-T regime. Then, the spin-phonon 
coupling will provide both the dissipation for phonons 
and the drag between phonons and spin excitations. 
In the drag conductivity, the phonon relaxation time 
(r ph oc l/|y b " ph | 2 ) enters together with the "drag rates" 
(l/ r ph-n-b K |^b-ph|2)_ As shown in Appendix 04, one 
can demonstrate that for quasiparticlcs with linear dis- 
persions and for k-independent boson relaxation time 
t£ = T b the following simplification for the drag con- 
ductivity is possible for an arbitrary form of the coupling 

T/b-ph . 
r p;k,k-p- 

P 

where u and v are the phonon and boson velocities, and 
Cph is the phonon specific heat. Since the diagonal con- 
ductivity of bosons in this case is K b = -p-v 2 r b Cb, where 
D s is the dimensionality of the spin system, the ratio of 
the drag conductivity to the boson one is independent of 
the scatterings and is defined by the boson and phonon 



specific heats 

Kb 3 Cb 

Since the population of phonons at a given temperature 
can be much larger than that of bosons, the drag conduc- 
tivity can significantly exceed the one by spin excitations. 
Similar argument is at the core of the original proposal by 
Gurevich for the large thermoelectric effect in metals,— ^ 
where the relation ftdrag/fc e = C p h/C e also implies the 
same drift velocities of phonons and electrons. 

While the parallel and the similarity between the 
electron-phonon drag and the thermal-only drag consid- 
ered in the last example are clear, they are not complete. 
The difference is in the presence of another diagonal con- 
ductivity term in our consideration, K ph oc l/|y b ~ ph | 2 , 
which, in the limit of the small spin-phonon coupling will 
dominate both Kdrag an d K b- Therefore, in general, the 
relation (|44|) does not imply the equivalence of the drift 
velocities of phonons and spin excitations. 



VII. CONCLUSION 

In this work we have considered a two-component sys- 
tem of phonons and spin excitations and have obtained 
general expression for the off-diagonal contribution to 
its thermal conductivity in the lowest order of the spin- 
phonon coupling. The off-diagonal contribution to the 
thermal current, referred to as thermal drag, is an en- 
hancement of the heat flux of one of the species due to the 
flow of another and vice-versa. We have employed two 
distinct approaches, the Boltzmann formalism and the 
Kubo approach, to derive the spin-phonon drag thermal 
conductivity and have established that both approaches 
yield identical results, Eqs. (fT7|) and (|39|) . In addition, 
we have considered contributions to drag from anoma- 
lous terms, which generally arise from the spin-phonon 
coupling, e.g. in the symmetry broken phases as well as 
in the gapped systems characterized by triplet-like ex- 
citations. While we mainly focus on the drag between 
phonons and bosonic spin excitations, we have also dis- 
cussed the case where the spin excitation's statistics is 
fermionic. 

To conclude, we have obtained an explicit expression 
for the drag conductivity in the two-component system of 
phonons and spin excitations under general assumptions 
on the nature of interaction between them. This should 
allow for the practical calculations of the drag effects in 
a number of materials. 
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Appendix A: Details of the Boltzmann approach 



(b) sC[n,f]=[ k > > , ] 



In this Appendix we provide some further details of the 
Boltzmann approach to the drag discussed in Sees. |TTT] FIG . 2: Graphical representation of the collision integral 
and El terms. 



1. Derivation of the drag rate of phonons on bosons 



Here we derive the "drag rate" of phonons on bosons 



ph->b 
p,k 



Eq. ©. 



While not necessary, it is nevertheless convenient to 
depict processes contributing to the collision integral as 
"probability diagrams" , see Fig. [5] The collision integral 
for bosons scattered off phonons via interaction ((3]) con- 
tains four terms, see Fig. [Ha): the first two increase the 
number of bosons with momentum k the other two scat- 
ter k-bosons into a different state. They can be grouped 
together by energy-conservation to yield, 



b-ph 
k-p I 



([/k- P np(/ k + 1) - /k(n p + l)(/ k -p + 1)](A1) 

x (5(e k - £k- P - w p ) 
+ [/k- P (n_ p + l)(/ k + 1) - / k n_p(/ k _ p + 1)] 
x (5(e k - e k _ p + Up) 

- P £- P ,k = (^k.k-p)* discussed after ©, writ- 
ing / = / + Sf and n = n° + Sn and neglecting terms of 
order SfSn and SfSf will yield the terms proportional to 
Sf and Sn shown in Eq. Performing this procedure 
and using <5n_ p = — Sn p yields the drag rate of phonons 
on bosons: 



using V 



Stt[f u ,n u + 5n] 



5n r 



/ j ph — >b : 
P T P,k 



ph-b 
r p,k 



= -2tt V 



b-ph 



p;k,k— p I 



(A2) 



[/k-p _ /k]^( £ k ~ £k- F 
-[/k°-p-/k°]^ k -^ 



Substituting this into the thermal conductivity in 
Eq. (flU)) and using relations ([To]) and ([To| yields the ther- 
mal conductivity in ([T4"]) and (|T7)) . 



The phonon collision integral discussed in Sec. IIIII B is 
shown in Fig. &h). 



Useful identities for the fermionic case 



Identities similar to Eqs. (fT5j) and ([16)) that are useful 
for simplifying expressions for the thermal conductivity 
and for relating its components to each other can also be 
obtained for the fermion-phonon system. They are 



dn° p 
duj 



p [tp-/^|[«M P + i] 



^ n P r fQ _ f0 



de 1 p 

5/k° 



(A3) 



w p-fk-^k- 



/k [ rO 

k-pj ~~ Q £ Uk-p 



(A4) 



^p — £k-p— £k 



These identities help to see that both contribution to the 
drag are equivalent. 



3. Derivation of the drag rates for the anomalous 
boson-phonon coupling 

The derivation of the drag in the case of the anoma- 
lous boson-phonon coupling is similar to the procedure 
detailed in Sec. IHII and Appendix [A"l 1 . For the coupling 
in (|41[) the scatterings describe the processes involving 
creation of two bosons from a phonon and generation of 
a phonon due to annihilation of two bosons. In that case 
the boson and the phonon collision integrals are described 
by two similar "probability diagrams" , see Fig. [31 

The expression for the boson collision integral has the 
following form, 

p 

x [(1 + /k- P )(l + /k)n p - /k- P /k(l + rip)] . (A5) 
Thus, the drag rate of phonons on bosons is: 
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St k [f,n]= 



ph r 

St p [n,f]= [ 



FIG. 3: Graphical representation of the collision integrals for 
the "anomalous" terms, Eq. (|41|l . 



s%[f 



br '° n° + 6n] « - 



y-r 5n p 



ph^b ' 
P T P,k 



1 _ Q_lf>b-Ph I 2 
ph-S-b ~~ Z7F I K P ;k,k- P | 



(A6) 



l + /£- P + /g)*(ek + ek-p-Wp). 



' P .k 



Similar consideration gives the phonon collision integral: 
Stf[nJ]=nJ2\V^^\ 2 (A7) 

ft 2 2 

k 

[(1 + ra P )/i=p/js±E - "pC 1 + /js^X 1 + /i±p)l 

22 2 2 

X (5(£ k- P + £ k+ P — OJp) , 

2 2 

where the factor of 2 has been removed to avoid double 
counting of the final states and the symmetrized nota- 
tions for the momenta are used. The drag rate of bosons 
on phonons is given by 



b— >ph ; 



1 



b->ph 



k T k,p 

^b-ph I 2 r „q 



'k,p 



x i5(ek + £k- P - w p ) 



(A8) 



After some algebra, the drag conductivity from bosons on 
phonons and phonons on bosons turn out to be identical 
and yield the total thermal conductivity of (|42|) . 



4. Drag conductivity in a limiting case 

Here we derive the thermal drag conductivity under 
two main assumptions: the boson scattering times are 
independent of the momentum and those of the phonons 
are determined entirely by its interaction with bosons. 
We consider linear energy spectrum for both bosons and 
phonons which are given by = u|k| and w p = u\p\, 
respectively. Thus the drag conductivity of Eq. (|17l) re- 



duces to, 



47ra 2 u 2 T b 



E k x x phlT/b-ph 1 2 
K P T p r P ;k,k-p| 



k,p 



S(s-k - £k- P - w p ) - <5(e k - £ k - P + w p ) 



(A9) 



The phonon scattering time in (|A9[) can be obtained 
from the phonon-boson collision integral in a standard 
way, similar to the derivation of the drag rates in Ap- 
pendix [A]l. Such a derivation yields, 



2 -Ei< k t P rx[/k-p-/k] (mo) 



8(ek - e k - P - w p ) - 5(e k - £k- P + w F 

We now rewrite the drag term in a compact form, 

Kdrag = -v 2 T h C ph) (All) 

using the auxiliary function Cph, which is a "modified" 
phonon specific heat given by, 



P 



(A12) 



The k-integration is now hidden in another auxiliary 
function F(p), which is given by 

vTF(p) = Ann? £V \V^f_ p \ 2 x [/°_ p - /°] 



S(e^ - £k- P - w p ) - <5(s k - £k- P + Up) 



(A13) 



Let us split the above expression into two terms 
k ~ I V pMM-p 

k 

xc5(e k -£k-p - Up), (A14) 



and 



xS 



^k- P 



(A15) 



In I2 we make the change k — p — >• k, 

/ 2 = -4^p h ^(k-- P -)|y p b ;p i_ k | 2 

k 

x[/£_ p -/£]d(e k -ek- P -Up), (A16) 
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Thus, C ph in (|A12|) is given by, 



J = dx/2-ir leads to 



p 



3u 2 ^ ^™P/„x\2 L^ph lT/ b -P h I 2 



T ^ doj 
p 



[j£-p-/ k Hek-Ek- P -< 



(A17) 



where we have utilized the relations \V, 



b-ph 



V 



b-ph 



2 and rj? 11 = r ph 



p;k,k— p I 

p k-pkl atxyjL 'p — '-p- Performing similar manip- 
ulations on tP' 1 in (|A10[) one can see that the k-integral 
in rP h cancels exactly the one in (|A17[) . Thus, 



P h 



3u 2 



p 



ph] 



(A18) 



where C p h is the phonon specific heat. This is a rather re- 
markable result given the arbitrary form of boson-phonon 
interaction. Thus, the thermal drag conductivity ex- 



pressed in terms of 



^drag 



; 2 r b 

Ds 
3 



C h /D s 

c b ' 



is 



(A19) 



where Cb is the boson specific heat, and D s is the dimen- 
sionality of the spin system. 

We find that Cph = C p h even for the scenario when 
3D-phonons interact with ID bosons. In this case, the 
3D-phonon momentum can be split into p = p|| +p_i_, the 
part parallel to the ID momentum of the boson k and 
the part perpendicular to it. Thus, the boson occupation 
numbers in the 3D-1D case change to (/k_ P|| ~fw) an( ^ * ne 

(5-function changes to S (ek~^k-p B + w p) ■ The momentum 
transformation to be used for this case is k — p|| — > k. 



n B (k,p,iw) = 4 



n(x)n(y) 



{-1m[G£(y)G2(x + y)] 

x G k (y - iu)G p {x + y - iu})G£_ p {-x) 
-Im[G£(y)GK(x + y)} 

x G k (y + ilo)G p (x + y + ilo)G£_ p (-x)} (Bl) 

+ {MG?(v)G2_ p (y-x)] 

x G k (y - ioj)G p {x - icj)G p (x) 
+G'i(y)G*{y + iu) 

x G k _ p (y - x + ilj)G p (x - iuj)G';(x)} (B2) 

+ {G£(y)G k (y-iu) 

x G k - P (y -x- iuj)G p (x + iuj)G p (x) 
+Im[G k i (y)G k i _ p (y - x)] 

x G k (y + iuj)G v (x + ioj)G';(x)}] 



n B (k,p^) + n^ c (k, P , 



<-) B <-) 

where G k W =Gk (iu. 



(B3) 
(B4) 



x + iO+) =g[ (x) + ih\ (x) 
refers to the retarded Greens functions and their decom- 
position into real (') and imaginary (") parts and n{. . .) 
is the Bosc distribution function. The subscripts 'a', 'b', 
and 'c' refer the contributions to n B (k, p, iui) which stem 
from the curly brackets labeled by Eqs. (|B1|) . (|B2[) . and 
(|B3|) . For the DC heat conductivity we need 



n B (k,p) = lim -lm\U B (k,p,iuj ^ uj + i0 + )] . (B5) 



We first take this limit focusing on Ii a B . The variables 
of integration can be substituted such as to express this 
limit in terms of a derivative of the distribution function 



Appendix B: Branch cut integration approach 



Here we derive the results of Section IIVI A using a 
different approach, namely by converting the frequency 
summations in (|2~5j) and (|2"r?|) to the problem of integra- 
tion along the branch cuts of the Green's functions. In 
Section IIVI A we have analyzed in detail the diagram 
A in Fig. [TJ Here we consider the derivation for the 
diagram B. Carrying out the Matsubara summation 
in Eq. (f2l)]) . obeying the location of the branch cuts of 
the Green's functions and using the shorthand notation 



n|(k, P ) 



n(x) 



dn(y) 



dy 

x{lm[G£(v)G%(x- 



u k- 



dn(y) ^ y,k,p) 
dy 



(B6) 



where the abbreviation i? a (x,y, k, p) has been defined. 
A similar substitution cannot be achieved for the contri- 
butions from Eqs. (|B2[) . and (|B3I) . Instead, we expand 
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the Greens functions to lowest order in u) 

d 



n b+c 



(k,p) 



n(x)n(y) 



x Jy 



duj 



GUy)G;(x) 

x {lm[Gf (y + oj)Gt P (y -x + u)G^(x - w)] 
-Im[Gf (y - w)Gf_ p ( y - x - u>)G*{x + u)]} 
+G';(x)Im[G£(y)G*_ p (y-x)] 
x{lra[G*(y + u)GK(x + u)} 

-MGk(v-u)G?(x-u>)]}]\ u = 



4/ / ra(x)n(2/)i? b+c (x,2/,k,p), 



(B7) 



As in Eqs. (|35|) . and ([36)1 we introduce a phcnomcnolog- 
ical, momentum dependent one-particle self-energy for 
the bosons [phonons] 



G[G](z,\l) 



1 



- £k[wk] + irk[sk]sgn(Im(z)) 



(B8) 



where [] -bracketed terms refer to phonons and z is com- 
plex. Inserting this into Eq. (|B6|) . R a (x, y, k, p) turns 
into a rational function 



R a (x,y, k,p) 



P a (x,y,k,p) 



(B9) 



Q Q (x,y,k,p) 
P a (x,j/,k,p) = -r k _p[r k w p (x + y) + s p £ k (y)] 2 
Q a (x, y, k, p) = [si + uj p {x + yff 

x [r k _ P + £k- P (-x) 2 ][7i + £ k (2/) 2 ] 2 , 

where the abbreviations w k (x) = x — w k and £ k (y) = y — 
£ k are used. Note that w k (x + y) = x + w k (y). Similarly, 
in Eq. (|B7[) Rb+ C (x, y, k, p) turns into 



-R&+c(x,y,k,p) 



P 6+C (x,y,k,p) 



(BIO) 



Qb+c(x,y,k,p) 
A +C (x,y,k,p) = 2s p |-2r k r k _ p [a; p (x) 2 + s 2 ] 
x [ r k + £k(y) 2 ]{r k -p[w p (x)r k - s p e k (y)] - [r k s p 
+w p (x)e k (y)]£ k _p(y - x)} + [r£_ p + £ k - p (y - x) 2 ] 
x (-(r k - r k _p) Sp [w p (x) 2 + s 2 ]£ k (y)[r 2 + e k (y) 2 ] 
+cj p (x)sp{-2r k r k - p + r k s p - r k r k _ p s p 
+r k s p e k (j/) 2 + 3r k r k _ p s p e k (y) 2 + 2r k _ p £ k (y) 4 
+4r k e k (y)[r k (r k + s p ) + £ k (y) 2 ]£ k - p (y - x)} 
+w p (x) 3 r k {r k - r k r k - P + 3r k _ p e k (2/) 2 
+?'k£k(y)[£k(?/) + 4£ k _ p (?/ - x)]})} 
Q b+c (x,2/,k,p) = [s 2 +w 2 (x)] 3 [r 2 +e 2 (y)] 3 

x [^k- P + 4- P (y-2;)] 2 



j-j-a(b+c) ^ ^ e evama t ec j assuming, as before, that 

the bosons and phonons are quasiparticles with r k (s k ) <C 
£k(wk)- In that case, expressions valid to leading order 
in r- k (s k )/e k (o; k ) for Eqs. (|B6|) . and (|B7|) are obtained 
from the residues of R a (b+c)(x, y, k, p) alone, while <zs- 
suming the distribution functions to be holomorphic and 
retaining only their lowest-order non-vanishing deriva- 
tives. Moreover, any imaginary part of the arguments 
of the distribution functions arising in that process 
can be dropped. Since the poles of R a (b+c) ( x , y, k, p) 
stem from quadratic equations at most, this calcula- 
tion can be done analytically. We emphasize that the 
proper evaluation of the higher-order contributions in 

^k(sk)/£k(^k) to Tl°^ b+C ^ (k, p) would require a treat- 
ment of the pole structure of the Bosc distribution func- 
tions and their derivatives. Analytically this is not fea- 
sible given R a (b+c) ( x i 2/? k, p). This also implies that the 
"non-Boltzmann" terms of Appendix [Clare not a system- 
atic account of all next-leading order corrections. The 
leading-order analytic calculation is tedious but straight- 
forward. After some algebra we arrive at 



n£(k,p) 



n^ +c (k, P : 



_L <M£k) 

2r k s p d£ k 

_L <9n(e k ) 

2r k s p <9£ k 



n(-s k - p )j-2 



?7k P 



n(w p ) 



k P c k P J 



k P 



kp 



(Bll) 



where 77 kp = r k + r k _ p + s p and e kp = uj p - e k + £ k _ p . 
Thus, the rightmost fraction in both expressions can be 
approximated by ir5(ui p — e k + £ k _ p ). The correspond- 
ing constraint lo p w e k — e k _ p has also been used to 
rearrange the arguments of the distribution functions in 
Eq. (|B11[) . Diagram A in Fig. [T]can be obtained directly 
from the preceding derivation by relabeling w 2 — > —1^2 
and by realizing that — iu n — lu p + is p sgn(hn(—iu n )) = 
— [iui n + ujp + is p sgn(lm(iui n ))]- I.e. xa can be ob- 
tained from Eq. (|B11[) simply by using the symmetries: 

' , and by replacing 
and the total drag is 



"A + K B 



^drag = + K21, the final result is 



- -— Wp, 1I X U, \V h ~ Ph \ 2 T b TP' 1 

g— rp 2^ W k £ kW p Wp| l/ p;k k _p| T k r p 



T 



k. P 



{^[l + /k°- P + <]^k-£k- 
-|^[/L p -n p ]^(£k-£k- 



- ^p) 



(B12) 



where we have renamed the Bose distribution functions 
with arguments £ k (w k ) to / k (n k ), as in section HVl A. 
This result is identical to Eq. (|39|) . Thus, both technical 
approaches within the Kubo formalism yield the same 
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Appendix C: Non-Boltzmann contributions 



In Section lTVl A we have discussed contributions from Im IT J 1 (oj) and showed that they lead to the results identical to 
the ones from Boltzmann theory. In the following wc will discuss additional contributions to drag thermal conductivity 

from the remaining terms of Eq. pip. hull™ (oj), m = 1,2 and 3, which, however, are subleading and can be safely 

neglected. Evaluation of these terms is rather cumbersome, and for illustrative purposes we focus only on Imll^ (oj), 
which is given by 



Imllj' (k,p,w) 



A k (e 1 )A k (e 2 )A k ^ p (e 3 )A p (e 5 +oj)A p (e 5 ) P- 



1 



£2 



(n 3 -ni)(n 3 _i -ns) (n 3 - ni)(n 3 -i - n 5+u ) , (n 3 - n 2 )(n 3 _ 2 - n 5 + u ) (n 3 - n 2 )(n 3 _ 2 - n 5 ) 



oj + e 5 + ei - e 3 cj + e 5 + e x - e 3 £ 5 + e 2 - £ 3 

where P stands for the principal value. In the limit of zero-frequency, oj — > 0, we obtain, 



£5 



£3 



Imn^ 2) (k,p,tj) 



j=0 



I I I I A k ( Sl )A k ( S2 )A k _ p (e 3 )(A p (e 5 )) 2 P-±- P-} 

J ei J e2 J E3 Je 5 V ' £l~£ 2 £1 



(ni - n 3 ) dn 5 

5 - £ 3 9£ 5 



(CI) 



TTT ■ (C2) 



Using the spectral representation (|35J) one can easily perform integrations in £ 2 . We further simplify the expression 
by performing the E\ and £ 3 integrations on the terms containing n 3 and n\, respectively, to obtain 



imni 2) (k, P ,w) 



j=0 



:(si)(A p (e 5 )y 



A 



A k - p (e 3 )[A p (s 5 ) 



1 



2/1 



9ni 

2y2T7f 9£7 

1 9 y 

2 dyy 2 + r k 



1 9 



2/1 



2 9yi yf + r£ 



9»5 
9£ 5 



9e s 



(C3) 



where y = (s k + £5 — £ 3 ) and yi = (si + £5 — £k- P )- The contributions from all three terms in the above expression 
are of the same order. Consider contributions from the first term which is given by 



Auie,) (A p (e 5 ) 



yi dni dn§ ir 1 



Vl 



dei de 5 



2 HcSp doj de 



ek-p). (C4) 



It appears that this expression contains the factor of the type x ■ S(x), which implies x = 0. However, under strict 
consideration, i.e., taking into account finite lifetime r^, x is non-zero and is of the same order as the "spread" of the 
5- function [x ~ r^). From a direct comparison of Eq. (|C4j) with Eq. (|37|) . we conclude that contributions from (|C4j) 
are smaller by the factor {r-k\df^/de\/f^) ~ rk/£k -C 1. Thus, the thermal conductivity contributions from (|C4|) and 
from the rest of the terms of (|C3[) can be neglected in comparison to the Boltzmann terms of Eq. (|37| . For the similar 
reason it is justified to use the delta-function form in Eq. p9[) for the leading contributions, because the broadening 
in the spectral function only yields a subleading correction of higher order in nc[s P ]. 
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